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The antiferromagnetic (AFM) order in the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) super- 
conducting state is analyzed on the basis of a Ginzburg-Landau theory. To examine the possible 
AFM-FFLO state in CcCoIns, we focus on the incommensurate AFM order characterized by 
the wave vector Q = Qo ± 9inc with Qo ~ (t, tt, tt) and qlnc || [110] or [110] in the tetragonal 
crystal structure. We formulate the two component Ginzburg-Landau theory and investigate 
the two degenerate incommensurate AFM order with glnc || [110] and [110]. We show that the 
pinning of AFM moment due to the FFLO nodal planes leads to multiple phases in magnetic 
fields along [100] or [010]. The phase diagrams for various coupling constants between the two 
order parameters are shown for the comparison with CeCoIus. Experimental results of the 
NMR and neutron scattering measurements are discussed. 
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1. Introduction 

The possible presence of a spatially modulated state in 
spin polarized superconductors was predicted by Fulde 
and Ferrell/' and by Larkin and Ovchinnikov^^ more 
than 40 years ago. While the Bardeen-Cooper-Schrieffer 
(BCS) theory assumes Cooper pairs with vanishing total 
momentum, the FFLO superconducting state represents 
a condensate of Cooper pairs with a finite total momen- 
tum. Since the FFLO state has an internal degree of free- 
dom arising from the reflection or inversion symmetry, a 
spontaneous breaking of the spatial symmetry occurs. 

Although this novel superconducting state with an ex- 
otic symmetry has been attracting much interest, the 
FFLO state has not been observed in superconductors for 
nearly 40 years. Naturally, the discovery of a new super- 
conducting phase in CeCoIns at high magnetic fields and 
low temperatures'^' triggered numerous theoretical and 
experimental studies because this high-field supercon- 
ducting (HFSC) phase is a most likely candidate for the 
FFLO state. The recent interest on the FFLO super- 
conductivity/superfiuidity extends further into various 
related fields, such as organic superconductors,^"^) cold 
atom gases,"'^"' ^"'^^ astrophysics, and nuclear physics. "'^^^ 

Although the HFSC phase of CeCoIns has been 
interpreted widely within the concept of the FFLO 
state,''' 13-21) recent observations of the magnetic order in 
the HFSC phase call for a reexamination of this conclu- 
sion.^^' Neutron scattering measurements have shown 
that the wave vector of the AFM order is incommensu- 
rate Q = Qo±9inc with Qo = (i, tf, tt) and the incommen- 
surate wave vector is fixed to ginc ~ (0.127r, ±0.127r, 0) 
independent of the orientation of magnetic field. ^'''^^' 
The AFM staggered moment Maf lit^s along the c- 
axis.^^'^^) For the magnetic field along [110] direction, 
the incommensurate wave vector is perpendicular to the 
magnetic field ginc ~ (0.127r, — 0.127r, 0), while the two 
degenerate incommensurate AFM states with q^^^ ^ 
(0.127r,0.127r,0) and (fine ^ (0.127r, -0.127r, 0) appear in 
the magnetic field along [100] direction. ^3' 



The magnetic order of CeCoIns appears in the super- 
conducting state, but not in the normal state. This cou- 
pling between the magnetism and superconductivity is 
in sharp contrast to the other heavy fermion supercon- 
ductors where the magnetic order is siipprcsscd by super- 
conduct ivity.'^^^ Another important feature of CeCoIns is 
the enhancement of HFSC phase under pressure.^'') This 
feature is also in contrast to conventional magnetic order 
in Ce-based heavy fermion systems which is suppressed 
by pressure. These unusual features indicate that the 
magnetic order of CeCoIns is not a conventional AFM 
order. 

Several intriguing quantum phases have been theoret- 
ically proposed for the unconventional magnetic order 
in CeCoIns.^^"'^^) Our proposal is based on the presence 
of AFM quantum critical point near the superconducting 
phase of CeCoIns.^^' '^^^ We have shown that the AFM or- 
der occurs when the inhomogeneous Larkin-Ovchinnikov 
state is stabilized in the vicinity of the AFM quantum 
critical point. ^^'^^^ The coupling between the magnetism 
and FFLO superconductivity seems to be consistent with 
the above-mentioned unusual feature that the magneti- 
cally ordered phase is confined in the superconducting 
phase in the H-T phase diagram. ^^'^"^^ Another proposal 
has been given on the basis of the emergence of a pair 
density wave state. ^^"^^^ 

In order to identify the HFSC phase of CeCoIng it 
is highly desirable to investigate these proposed phases 
for a comparison with experimental results. For this pur- 
pose, we investigate the AFM phases in the FFLO state. 
We show that the multiple phases appear in the AFM- 
FFLO state when the magnetic field is applied along [100] 
or [010] direction. We discuss the possible phase diagram 
of CeCoIns on the basis of the comparison with experi- 
ments. 

2. Formulation 

An intriguing magnetic phase diagram arises from the 
degeneracy between the incommensurate wave vector 
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(fine (0.127r,0.127r,0) and gi„c ~ (-0.127r, 0.127r, 0) 
in the magnetic field along [100] or [010] direction. In 
order to investigate the magnetic phases, we formulate 
the phenomenological two component Ginzburg-Landau 
model. We here consider the AFM order in the inhomoge- 
neous Larkin-Ovchinnikov state for which the Ginzburg- 
Landau functional of the free energy is described as, 

F{vi, V2)/Fo = [{T/T° - 1) + eAFiQi - &]vl 

-\'nim^C2{n)5{qi - q2,2nqFFLo), (1) 

n 

where rji and 772 arc the two component order parameters 
corresponding to the two degenerate AFM states with 

Q = Qo ± (fi ~ Qo ± qinJ and Q = Qo±q2^ Qo± Qi^J, 
respectively. The wave vector of incommensurate AFM 
order is assumed as q^^J = (0.1257r, 0.1257r, 0) and q^^ = 
(— 0.1257r, 0.1257r, 0). Small deviations of wave vectors 
from q^^J and q^l^J due to the pinning of FFLO nodal 
planes are taken into account in eq.(l). In our study the 
incommensurate wave vectors q^^J and q^^^J are not mi- 
croscopically derived but assumed on the basis of the 
experimental results in Refs. 23 and 24. We think that 
the q^^J and q^l^J are pinned through nesting features in 
the band structures. 

We define the AFM staggered moment MApir) = 
(— l)^+^+^M(r) where M{r) is the magnetic moment 
perpendicular to the applied magnetic field at r = 
{x,y,z). Then, the AFM staggered moment is described 
as, 

MApir) = Mo[?7i cos(gi ■r)+r]2 cos(g2 • r)], (2) 

where Mq (and Fq) are scaling factors by which the 
Ginzburg-Landau free energy is written in the renormal- 
ized form (eq.(l)). We take the unit of magnetic moment 
and energy so that Mq = Fq = 1. The is the Neel 
temperature for ci = C2 (n) = and ^af is the coherence 
length of AFM order. 

The phase diagram of the Ginzburg-Landau model 
considerably depends on the sign of the non-linear quar- 
tic coupling term b. The single-g magnetic structure 
{111,112) oc (1, 0) or (0, 1) is stabilized by the positive cou- 
pling constant b, while the negative b favors the double-g 
magnetic structure (r?i,??2) cx; (1,±1). 

The coupling constant ci describes the effect of mag- 
netic field H = (_ffx, Hy, Hz) which may break the degen- 
eracy of 771 and ri2- According to the neutron scattering 
measurement for H \\ [110],^^^ ci is positive in GeCoIns. 
This is consistent with our theoretical analysis based 
on the Bogohubov-de-Gennes (BdG) equations.^^'^'^'^^^ 
When we consider the magnetic field along [100] direc- 
tion, this term vanishes and therefore the degeneracy of 
rji and 772 remains. We focus on this case in this paper. 

Effects of the broken translation symmetry in the inho- 
mogcneous Larkin-Ovchinnikov state are taken into ac- 
count in the commensurate term, that is the last term of 



eq.(l). We define 6{qi — 92, 2n(fFFLo) = 1 when the com- 
mensurate condition q\ — q2 = 2ngFFLO is satisfied for 
an integer n, and otherwise 5{qx — q2, 2ngFFLo) = 0. The 
modulation vector of FFLO state is denoted as ^fflO) 
and then the order parameter of superconductivity is de- 
scribed as A(r) — Aosin((fFFLO • r)- The commensurate 
term describes the pinning effect of FFLO nodal planes 
for the AFM moment. According to our analysis based on 
the BdG equation, the AFM moment is enhanced around 
the FFLO nodal planes where the supc;r(;onducting order 
parameter vanishes.^^'^'''^^) Then, we obtain the positive 
coupling constant C2{n) > 0. 

3. Magnetic Phases 
We here consider the FFLO state in which the higher 

harmonic FFLO wave vector 2Nqppi,o is close to q^^J — 

(2) 

^inc • Then, the broken translation symmetry plays an 
important role for the appearance of magnetic phases 
as will be shown below. When the wave vector q^^^J — 

q^l^ is not in the vicinity of 2n(fFFLO for any integer 
n, the commensurate term can be neglected. Then, the 
singlc-g magnetic phase with (771,772) oc (1,0) or (0, 1) is 
realized for a positive 6, while the double-^/ phase with 
(771,772) oc (1, 1) is stabilized for a negative b. We choose 
the parameter C2{N) — 0.01 in the following part and 
investigate the magnetic phase diagram for various b. 

We determine the magnetic phase by minimizing the 
Ginzburg-Landau free energy (eq.(l)) with respect to 
the order parameters 771 and 772 and their momentum 
gi and 5*2 for each temperature T and FFLO wave vec- 
tor qpFLO- We assume 771 > 772 > without any loss 
of the generality. As for the direction of (fpFLO) we as- 
sume qpFLO = qpFLOx for H \\ [100] as in Refs. 27 
and 35. Then, the phase diagram is determined by the 
renormalized parameters T/T^ and ^af^Oj where the pa- 
rameter qo describes the mismatch of the higher har- 
monic FFLO wave vector 27VgFFLO and the incommen- 
surability along the i-axis q^^^J — q^l^ ■ We define go as 

QinJ - ^inJ + 2^0^ = (2ginc + 2go)^ = 27VgFFLO with 
(Zinc = 0.1257r. Note that go = when the commensu- 
rate condition q^^J — q^l^ = 2AfgFFLO is satisfied. Since 
the FFLO wave number gFFLO increases with increasing 
the magnetic field and/or decreasing the temperature, 
the parameter go increases with the magnetic field, and 
change its sign at the commensurate line in the H-T 
phase diagram on which the commensurate condition is 
satisfied. 

We find that three magnetic phases are stabilized. One 
is the single-g phase where 

Vi>V2, (3) 

9*1 = + (1 - x)qoX, (4) 

Q2 = Qi^J - (1 + x)qoX. (5) 

The mirror symmetry with respect to the x- and y-axes 
is spontaneously broken in the single-g phase. The others 
are the double-g phase and double-g' phase where 

Vi=V2, (6) 
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Fig. 2. (Color online) (a) Phase diagram and (b) the ratio of order 
parameters 7?i/r)2 for 6 = 0.01. 
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Fig. 1. (Color online) (a) Phase diagram of the Ginzburg-Landau 
model with b = 0.1 for ^af<?0 a-nd the renormalized temperature 
T/T°. The definition of qo is given in the text. The go increases 
with the magnetic field, (b) The ratio of order parameters ri\/ri2. 



The analysis of the quadratic terms in eq.(l) shows 
that the double-g phase is stabilized immediately below 
the Neel temperature for (^af9o)'^ < C2(iV)/8, while a 
single-g phase is stabilized for (^af9o)^ > C2{N)/8. The 
Neel temperature is increased around the commensurate 
line ^af9o = owing to the broken translation symmetry. 
However, the enhancement of the Neel temperature Tn ~ 
J.0 _ C2(^jv) jg negligible unless the coupling constant 
C2{N) is large. 

Figures l(a)-4(a) show the phase diagram obtained by 
the numerical calculation for various coupling constants 
6, while Figs. l(b)-4(b) show the ratio r]i/i]2- First, we 
discuss the magnetic phases for a positive b. We see that 
the single-g phase is stable at low temperatures in Figs. 1 
and 2. This is because the quartic term 677^77! favors the 
single-g phase. Figure 1 shows that the phase diagram is 
mostly covered by the single-g phase for a large coupling 
constant b » C2(iV). Then, the double-g phase appears 
around the Neel temperature when I^af^oI ^ \/c2(iV)/8, 
but it disappears by decreasing the temperature. The 
double-g phase is stabilized in the intermediate tem- 
perature region T ^ O.STn for a small and positive 



b = C2{N) (Fig. 2). Figures 1(b) and 2(b) show the ra- 
tio r]2/r]i decreases in the single-g phase with decreas- 
ing temperature T/T^ and/or increasing the mismatch 
|^af9o|- The phase transition between the single-g and 
double-g phases is second order. 

Next, we discuss the magnetic phases for a negative 
coupling constant b. Figures 3(a) and 4(a) show that 
the double-g phase is stable at low temperatures for 
ICaf9o| < \/c2(iV)/4 while the double-g' phase is stabi- 
lized when the mismatch is larger than the critical value 
ICaf9o| > \/c2iN)/4:. The single-g phase appears around 
the Neel temperature near |^af9oI = \/c2(^)/4. It is 
shown that the phase diagram is mostly covered by the 
double-g and double-g' phases. The phase transition is 
first order between the double-g' phase and the other 
phases, while that from the double-g phase to the single- 
q phase is still second order. 

In order to understand the magnetic phase diagram 
more clearly, we plot the <^af<?o dependence of order 
parameters for b ~ 0.1 in Figs. 5(a)-5(c) and that for 
b = —0.01 in Figs. 5(d)-5(f). It is shown that the mag- 
netic phases are almost independent of b around the 
Neel temperature (Figs. 5(c) and 5(f)), while the single-^ 
phase (double-g or double-g' phase) is stable at low tem- 
peratures for the positive (negative) coupling constant b. 
Fig. 5(e) shows a discontinuous jump of order parame- 
ters at the first order transition from the single-g phase 
to the double-g' phase. 

4. High Field Phase Diagram of CeCoIns 

We discuss the possible AFM-FFLO state of CeCoIns 
on the basis of the Ginzburg-Landau theory given in this 
paper. For this purpose, we show the schematic phase 
diagram of CeCoIns for the magnetic field and tempera- 
ture in Fig. 6. 

The maximum amplitude of the FFLO modulation 
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Fig. 5. (Color online) ^af90 dependences of the order parameters 
rii and »?2 for b = 0.1 (Figs. 5(a)-5(c)) and b = -0.01 (Figs. 5(d)- 
5(f)). T = 0.5T° in Figs. 5(a) and 5(d), T = 0.9T° in Figs. 5(b) 
and 5(e), and T = 0.99T0 in Figs. 5(c) and 5(f), respectively. 



Fig. 3. (Color online) (a) Phase diagram and (b) the ratio of order 
parameters ri\/r)2 for b = —0.01. The dashed lines in (a) show 
the first order phase transition. 
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Fig. 4. (Color online) (a) Phase diagram and (b) the ratio of order 
parameters r;i/??2 for b = —0.1. The dashed lines in (a) show the 
first order phase transition. 



vector (7FFLO is approximately the inverse coherence 
length (/FFLO ^ According to the experimental es- 
timation of the coherence length of the superconductiv- 
ity ^,1^) the minimum number of N which satisfies the 



commensurate condition q-^^^ — q^^-^^J = 2Af(fFFLO is ap- 
proximately obtained as ^ 4. Then, the commensu- 
rate condition is satisfied in the FFLO state on a se- 
quence of commensurate lines (dashed lines in Fig. 6) 
with TV = 4, 5, 6, 7 



The magnetic phase diagram in the FFLO state is 
quite different between the positive and negative cou- 
pling constants b, as shown in Figs. 6(a) and 6(b). Fig- 
ure 6(a) shows the phase diagram for a positive b. We see 
that the AFM-FFLO state is mostly the single-g phase 
while the double-g phase is stabilized around the AFM 
transition line near the commensurate lines (shaded area 
of Fig. 6(a)). Note that the double-g phase does not ap- 
pear around the first order normal-to-FFLO transition 
line at which the AFM moment as well as the super- 
conducting order parameter are discontinuous. Since it 
is expected that the coupling constant C2{N) decreases 
with increasing N, the double-q phase is suppressed with 
decreasing the magnetic field. 

Figure 6(b) shows a schematic phase diagram for a 
negative coupling constant 6, where the double-q and 
double-g' phases are stable at low temperatures. In con- 
trast to Fig. 6(a), a sequence of discontinuous first order 
transition lines appear in the AFM-FFLO state. These 
first order transitions occur due to the pinning of AFM 
moment on the FFLO nodal planes. This is a direct con- 
sequence of the broken translation symmetry in the in- 
homogeneous Larkin-Ovchinnikov state. However, these 
first order transitions have not been observed in the ex- 
periments of CeCoIns. This indicates that the coupling 
constant b is positive in CeCoIns. Thus, the phase dia- 
gram in Fig. 6(a) is more likely realized in CcCoIns than 
Fig. 6(b). Since the second order phase transition shown 
in Fig. 2 has not been observed in the magnetic phase of 
CeCoIns, The coupling constant b should be much larger 
than C2{n). Then, the AFM-FFLO state is mostly cov- 
ered by the single-g phase where the mirror symmetry 
with respect to the x- and y-axis is broken. It is desirable 
that the phase diagram of CeCoIns will be investigated 
by the intensive experiments in the magnetic field along 
[100] axis. 

The phase diagram in Fig. 6(a) is supported by the 
nuclear magnetic resonance (NMR) and neutron scatter- 
ing measurements for CeCoIns. We first discuss the neu- 
tron scattering measurement. According to the results 
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(a) b>0 



AFM-FFLO 




(b) b<0 

HI AFM-FFLO 




Fig. 6. (Color online) Schematic phase diagram of CeCoIns for 
the magnetic field along [100] direction, (a) 6 > and (b) fe < 0. 
"BCS", "FFLO", and "AFM-FFLO" states are shown in the 
figure. The dashed lines show the commensurate lines where the 
commensurate condition is satisfied. The shaded region shows 
the double-q phase in (a), while it shows the single-g phase in (b). 
The dash-dotted lines in (b) show the first order phase transition 
lines between the double-g phase and double-g' phase. 



obtained by Kenzelmann et al.,, the position of Bragg 
peaks does not change with increasing the magnetic field 
along [too] direction.^'') The phase diagram for a nega- 
tive h (Figs. 3 and 4) is incompatible with this experi- 
mental result. Figure 7 shows the ^af9o dependence of 
the shift of main Bragg peak Agx = <fi — q-^^^ ■ For 6 < 0, 
Fig. 7(b) shows the shift ^ApAq > 0.05 at low tempera- 
tures. When we assume ^af = 3, Aq ~ O.OOGtt near the 
first order transition lines. This shift is larger than the 
experimental error, ^^-^ and therefore can be observed if 
it would occur. However, the experimental result has not 
shown the shift of Bragg peaks. 

On the other hand, the neutron scattering measure- 
ment is consistent with our results for the positive h. 
Fig. 7(a) shows a pronounced shift ^af Aq up to 0.03 near 
the Neel temperature, but the shift rapidly decreases 
with decreasing the temperature. Since the magnetic mo- 
ment is tiny near the AFM transition, it is difficult to ex- 
perimentally observe of the shift of Bragg peaks around 
Tn. Thus, our results in Fig. 7(a) are consistent with the 
experimental observation.^^' 

We here comment on the effect of domain structure in 
the single-g phase. For the magnetic field along [100] di- 




(b) b=-0.01 




Fig. 7. (Color online) The shift of the position of main Bragg 
peak from gt^'. We plot CafAi} for T = O-ST^, T = 0.9T§,, and 
T = O.ggrO, where Aq is defined by Aqx 
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rection the elastic Bragg peaks appear at Q — Qq ± q^^J 

as well as at its symmetric point Q = Qq ± 'finc^-^'^'^^^ 
This result seems to be incompatible with the phase di- 
agram for 6 > 0, because the intensity of Bragg peaks at 
Q = Qo i q^^^ is very weak in the single-g phase with 
Ivil ^ l'?2|- However, the observed four Bragg peaks are 
consistent with the single-g phase by taking into account 
the domain formation of two degenerate single-q states, 
i.e. \t]i\ > I772I and |?7i| < |7y2|- 

The phase diagram in Fig. 6(a) is also supported by the 
NMR measurements which can distinguish the single-q 
phase from the double-q phase. When the hyperfine cou- 
pling has the dipolar symmetry, the NMR spectrum at 
the In(2b) site shows the distribution function of the in- 
ternal field arising from the AFM staggered moment. 
The experimental results show the double peak structure 
of NMR spectrum at the In(2b) sitc^^'^^^ 37-39) Accord- 
ing to our analysis which will be published elsewhere, 
these results are consistent with the single-g phase, but 
not with the double-g phase. Thus, both neutron scat- 
tering and NMR measurements support the positive cou- 
pling constant b (Figs. 1 and 6(a)) and rule out the neg- 
ative b (Figs. 3, 4 and 6(b)). 

When we take into account the broken translation 
symmetry arising from the vortex lattice, other com- 
mensurate lines can appear in the phase diagram and 
stabilize the double-g phase around T — T^. However, it 
is expected that the effect of vortex lattice on the mag- 
netic order is smaller than that of the FFLO nodal planes 
when the Maki parameter is large enough to make the 
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lattice spacing of vortices much larger than the coherence 
length. 

5. Summary 

We investigated the incommensurate AFM order in the 
FFLO superconducting state with a particular interest 
on the possible coexistence of the AFM order and FFLO 
superconductivity in the HFSC phase of CeCoIns. As- 
suming the incommensurate AFM order with (fine || [110] 
or 5inc II [110]; we examined the magnetic phase in the 
magnetic field along [100] direction. We find that the 
multiple AFM phases appear in the H-T phase dia- 
gram owing to the broken translation symmetry in the 
FFLO superconducting state. The comparison between 
the Ginzburg-Landau theory and experimental results 
shows that the AFM-FFLO state is mostly covered by 
the single-g phase in which the mirror symmetry with 
respect to the x and y axes is broken. 
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